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Abstract. We study the elasticity of the plane-filling Peano chain in terms of a spring 
network, modeling a flat and thin elastic strip, smoothly curved so as to follow the 
chain pattern. 

This network sustains normal modes, named Peano modes; in the same way as acoustic 
waves in crystals are sustained by the matching of translations with rotations, these 
modes arise when discrete dilatations match with 7r/4 rotations. The pattern shares 
this eightfold symmetry with the octagonal quasicrystal, but is much looser and prone 
to disruption; its relaxation marks the saddle point separating the solid from the liquid. 
In terms of the equivalent quantum mechanical problem the above corresponds to the 
D=2 quantum delocalization transition. 

Some experimental observations, regarding a priori uncorrelated systems such as 
colloids, bilayer water and chromosomal DNA come together, at least qualitatively, 
within the viewpoint proposed here. With reference to the D=2 colloids the model 
accounts for the bosonic peak found in their phononic spectrum and for the relaxation 
of the excitations. The latter typically evolve towards patterns of increasingly large 
size, as in the inverse energy cascade of D=2 turbulence. As for bilayer water, it 
undergoes a liquid to quasicrystal transition surprisingly similar to the phase transition 
found here. Finally, considering the chromosomal DNA, the model unveils the mixed 
nature of its compact conformations, lying at the boundary between a liquid and a 
solid, similar to colloids and quasicrystals. 
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1. Introduction 

In this work we study the harmonic spectrum and the time relaxation of the planar 
Peano-Hilbert curve PH 2 . The curves PHo that £11 D-dimensional spaces are used in 
ordering multidimensional data [1, 2], as well as in manufacturing miniature antennas 
[3]. In the domain of polymer physics the curve PH^ has been proposed as a possible 
final state of a scale-invariant collapse [4], Quite remarkably, given the high degree 
of ordering of PH^, its pattern has been found to fit with the spatial organization of 
the chromosomal DNA in humans over a rather large range of scales [5]. This resulted 
from detecting the physical contact frequencies between pairs of genomic loci, using 
chromosome conformation capture methods [6]. Regarding DNA dynamics, significant 
data emerged from experiments on bacterial (E.coli) chromosomes [7, 8]. It was found 
that DNA loci undergo anomalous diffusion with a rather sharp distribution in the 
exponents, where the dominant behavior can be connected with a fractal dimension 
Df < 3 [9]. A physical description of the chromosomal conformation could be of interest 
for cell biologists and should be based on quantitative informations regarding the chain 
topology, the distribution of crosslinks, the nature of interchain interactions; in addition, 
this looked-for description should account for the complex multicomponent medium of 
the cell. Originally motivated by this phenomenology, the present paper is limited to 
analizing the harmonic properties of the planar curve PH 2 . We consider a quasi one¬ 
dimensional elastic strip, modeled as a mechanical network running along PH 2 , where 
deviations from the exact pattern carry a harmonic energy. The elastic couplings are 
anisotropic since they vary according with the local geometry, so that the resulting 
mechanical network turns out to be rigid: apart from uniform translations it does not 
allow for zero-energy modes. In particular, rigid rotations or uniform dilatations of any 
portion of the network have an energetic cost. The model does not account for the 
interaction among physically adjacent sites when they are distant along the chain: our 
aim is to describe waves propagating along meandering one-dimensional patterns, not 
necessarily having the nature of covalent bonds. This description could as well apply to 
electromagnetic waves, as it is known that the D=2 elasticity has strong similarities with 
the D=2 electromagnetism: the elastic deformations can be decomposed into irrotational 
and solenoidal components, where dilatation (D) and shear {S) have the roles of electric 
and magnetic field respectively. The main difference between the two theories regards 
the wave propagation properties, because dilatation pays an higher energy with respect 
to shear, and as a consequence in the bulk the longitudinal (D) field propagates with an 
higher velocity with respect to the transverse {S) field. The two distinct waves couple at 
the medium’s boundary, where they form the Rayleigh waves, which are characterized 
by nontrivial boundary-dependent dispersion laws. This dependence is exploited to 
extract information during earthquakes, and in the case of electromagnetism it can 
be used both to extract and to manipulate information, as done, e.g., by scattering 
electromagnetic waves through photonic crystals [10] or by engineering e.m. wavefronts 
with metamaterials [11]. In the elastic model in question here the particular ordering 
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enforced by the Peano pattern, which is built by inflating quadrupolar partitions of 
planar regions, produces non-dispersive Rayleigh waves, named Peano waves. 

In order to make clear the physical interpretation of the network, we recall how contact 
interactions among rigid disks undergo harmonic modeling, a procedure going under the 
name of fluctuation method [12], In that framework the notion of contact (harmonic) 
network [13, 14] can be given the following precise meaning [15, 16]: ’In a metastable 
state of a hard sphere system, two particles are considered to be in contact when they 
collide during some time interval At much smaller than the relaxation time of the 
structure, where metastability is lost, and much larger than the (thermal) collision time’. 
If the Peano pattern is to be associated with a linear polymer, the model describes the 
polymer’s fluctuations within a steep effective potential. Notice that the potential has its 
own time evolution, but on the larger scale of the structural relaxation. At shorter times 
the line of minima of the potential follows the topology assigned by the ordered sequence 
of bendings of the undeformed pattern. Within this polymeric interpretation, if one were 
to add to the picture the contact couplings among adjacent but nonsubsequent sites, one 
would end up with an elastic network describing a textured membrane. The phononic 
spectrum of such enlarged model would then display in the density of states, in addition 
to the contributions described in this work, the behavior p{\)d =2 ~ A, as expected 
from standard D=2 Debye spectra. It is worth mentioning that some experiments on 
two-dimensional colloids did in fact display an anomalously high density of states at low 
frequencies [17, 18], superimposed to the D=2 Debye behavior. As a matter of fact we 
show here that the connection is not limited to the so-called bosonic peak of colloids, 
but involves their space-time behavior as well. In analogy with quasicrystals, whose 
phononic spectra are marked by self-similarity [19, 20, 21], we exploit the recursive 
nature of the curve and consider the sequence of its periodic approximants having the 
form of closed chains on the square lattice. At the lattice sites the harmonic operator acts 
on the two-dimensional space of planar deformations; it has a (2 x 2) block tridiagonal 
structure since it involves, at each site, the two nearest neigbours along the chain. For 
the sake of brevity we will call this operator disordered laplacian, in spite of it being 
neither scalar nor disordered. In fact the operator is deterministic, as the sequence of 
the different vertices that assigns the chain pattern. 

We have found that in the case of PH 2 the spectral problem for the disordered laplacian 
translates into a hnite-dimensional trace map, which can be treated by standard 
methods. The phononic spectrum of the closed curve shows a multiplicity of power- 
law singularities in the density of states, with exponents covering a whole interval (see 
Figs.l and 2). The physics implied by such an exotic spectrum emerges from the analysis 
of the relaxation process, described by the overdamped Langevin equation associated 
with the disordered laplacian. The result is that every exponent identifies the scaling 
behavior (with respect to length and time) of a specihc class of fluctuations. 

We then found that the anisotropy of the elastic coupling induces a phase separation 
(see Fig.3); the two phases are respectively dominated by dilatation {V) and shear 
(iS), and can be identihed with solid and liquid. The critical point where the two 
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phases coalesce corresponds to the isotropic limit; quite remarkably the interference 
of (V) and (S) generates not only anomalous fluctuations, but also phononic modes, 
the Peano modes. During the relaxation process the energy flows from localized to 
extended structures, as in the D=2 turbulence [22] and as observed in the mentioned 
experiments on two-dimensional colloids [17, 18]. From the perspective of the quantum- 
mechanics determined by the disordered laplacian, the isotropic limit corresponds to 
the critical point of the D=2 quantum delocalization transition. In other words the 
Peano chain identihes a two-dimensional extension of the D=1 Aubry-Andre’ class [23]. 
A comprehensive exploration of the critical point requires an analysis of the trace map; 
here we limit ourselves to some preliminaries regarding its symmetry properties. As 
discussed in Section 6 we hnd an eightfold symmetry, revealing that PII 2 is close to 
the octagonal quasicrystal, but with a definitely lower connectivity. Interestingly, a 
liquid to dodecagonal quasicrystal transition has been reported for bilayer water [24] 
in a similar regime. Some preliminaries regarding the properties of the Peano chain 
are reported in the Appendix where, more importantly, we derive a substitution rule for 
closed chains. In Section 2 we introduce the network model and the disordered laplacian; 
we then formulate the eigenvalue problem for this laplacian in terms of a trace map. In 
Section 3 we analize the singularities of the spectral density of states by means of the so- 
called thermodynamic formalism of multifractals. Readers uninterested in the spectral 
properties can skip Sections 2 and 3 and go over to Section 4, where we examine the 
fluctuations. The anomalous diffusion exponents and fractal dimensions are explicitly 
obtained for the planar case and estimated by means of a mean-held argument valid for 
chains of arbitrary dimension D. The exponents obtained with D = 3 appear to be in 
agreement with measures of anomalous diffusion performed on bacterial nucleoids [7, 8]. 
Section 5 is centered on the solid-liquid phase transition and on the properties of the 
relaxation process. In Section 6 we discuss the symmetries of the system at the critical 
point and illustrate how they relate with the Peano modes. 

2. Elastic network and trace map 

We consider the vector held of deformations of the plane-hlling Peano-Hilbert chain 
PH 2 . An example of the patterns described by PP [2 is given in Fig.2A; as it is shown in 
the Appendix, the closed chain can be constructed by means of an iterative procedure 
dehned over the square lattice. One starts with the approximant of order k = 1, taken 
as the square path, and inhates the paths with suitable rules. The approximant of order 
/c is a closed chain that covers with no intersection all sites of the (2^ x 2^), [k = 1, 2,...) 
lattice. Upon assigning an orientation to the path the lattice sites can be labeled with 
the chain coordinate n, {n = 1,2, ...4^), and the vector held of deformations can be 
written as x{n) = {x{n),y{n)). This labeling is appropriate in our context since we 
assume that the excitations are transmitted only along the Peano pattern and that 
their elastic energy depends on its local geometry. The dynamics studied here involves 
the two-dimensional spaces of deformations at the sites (n — 1, n, n -|- 1). The normal 
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modes of the model solve the following equation: 

A ix(n) = - 1) + + 1 ) + 14^(n), (1) 

where I is the identity and ^(out) 2 x 2 matrices describing the contributions 

from the (incoming/ outgoing) links attached to the n-th site. The general form of the 
T’s, when including an interaction that mixes the x and y components, is given by: 

^x,x kxj ^x,y q (2) 

'^y,x = 7 ) '^y,y ~ 

The operator V is the contribution from the vertex that is determined by requiring the 
invariance under uniform translations: 

14 + fi'”) + = 0, (3) 

this condition thus guarantees that the eigenvalue A = 0 is in the spectrum at its upper 
bound. The network is anisotropic because we assign an higher energy to the deviations 
when they are orthogonal to the local alignement of the pattern. If we label with a t 
(transverse) and with an / (longitudinal) the corresponding elastic constants we thus 
have kt, ki, kt > ki. In correspondence with an horizontal link the diagonal entries 
of the matrix T are k^ = h, ky = kt] the off-diagonal entries 7 can describe, in the 
polymeric interpretation, the elasticity of the chain bendings. One can notice that even 
in the absence of this contribution the model has no soft (zero energy) modes other 
than the uniform one. On the r.h.s. of the eigenvalue equation one has an operator 
with the appearance of a matrix-valued inhomogeneous laplacian that acts on the two- 
dimensional vector held of deviations. For the sake of brevity let us name it disordered 
laplacian A(^ph), with the proviso that its inhomogeneity is actually deterministic since 
it is determined by the Peano pattern. When the off-diagonal contribution is disregarded 
(7 = 0 ) the X and y components are decoupled, and one is left with two scalar disordered 
laplacians A^, and Ay-. 

A(ph) = Diag{Ax{kt, h), Ay{kt, h)). (4) 

These operators are essentially identical because they turn one into the other under 
exchange of the couplings {kt ki, ki ^ kt): 

Ay{kt,ki) = Ax{ki,kt), (5) 

and for this reason we will focus our analysis on A^,. Let us thus write the equation for 
the component x{n,t): 

(A - Vn)x{n) = T^t'^x{n - 1) + T^r^^x{n + 1) ( 6 ) 

(-p{out) _ pC*"-) \ 

I-*- n ~ n+l) 

Here the coefficients T are c-numeric and have value ki{kt) at the horizontal (vertical) 
links; correspondingly the site (vertex) potential is given by 14 = By 

following a standard procedure, this equation can be written in the form of a single step 
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map acting on a two-dimensional space. Once the transfer matrix T = T(n) has been 
assigned to the n-th vertex, the vector z{n) = {x{n),x{n — !))"'' is mapped by T into 
z{n + 1): z{n -|- 1) = T{n)z{n). The transfer matrix T = T{n) has the form; 




(A - Vn) 

■p{out) ’ 

i n 


T { n ),,2 = 


■p(m) 

n 

■p{out) 

i n 


( 7 ) 


T(n)2,i = 1, T(n)2,2 = 0 


In this way the eigenvalue problem translates into a boundedness condition for the 
transfer matrix Ttot of the total chain Ttot = T{Nk) ■ T{Nk — 1) • .... ■ T(2) ■ T(l). It 
is convenient to make use of SL{2, R) (i.e. determinant -|-1) matrices, because the 
boundedness condition for Ttot can then be written as a single condition for the trace: 
l\Trace{Ttot)\ < 1- 

In view of that we perform the following transformation: 

T(n) ^ f(n) = + 1) ■ T(n) • ^(n), 5(n) = 

In fact the new matrices T = Tout,in = T{a, (3) belong to SL{2, R) and their entries are: 

Ti 1 = "x + H—; Ti 2 = —O' (8) 

p a 

^ 2,1 = 1 /^ 2,2 = 0 


The parameters a, /3 introduced above are functions of the hopping coefficients: a = 
and j3 = /9(n) = The iteration formula obtained 

in the Appendix (see Eq.A3), that we are going to represent by means of transfer 
matrices, involves an operation T that describes the inversion of the ordering along the 
chain. When applied to a single step transfer matrix this inversion has the effect of 
exchanging -w-ith T^*”) because obviously it must be T(To«t,m) = Tin,ouu in terms 

of the parameters a, /3 we have: 

r(f (a, /3)) = f (1/a, /3) = ai ■ f-\a, /3) ■ a,, (9) 

where Ui is the Pauli matrix [(o'i)i 2 = (o'i) 2 i = 1, (c^i)!! = (o'i )22 = 0]. The action of T 
on multiple step transfer matrices M G SL{2, R) is thus given by T(M) = ai ■ 

We label the transfer matrices of the vertices with the directions of their links as follows: 
Txy, Txx, Tyx, Tyy, whoro the left (right) label refers to the outgoing (ingoing) link. The 
inflation rule for Tyx and Tyy has the form: 

f'yX = Tyy ' Tyx ' T {T^y ' Tyx) (lO) 

f'yy=fyx-fxyT{fyx-fxy), 

the ’conjugate’ vertices T^y and T^x transform according to the analog of this formula 
where the labels x and y are exchanged. 

By using the representation of T we have: 

f'xy = Tyy ■ fyx ' ' T ' f-^ ' 


( 11 ) 
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'^'yy ~ ' '^^y ' *^1 ' '^xy ' '^yx ' *^1 • 

From this equation we computed the band spectra on closed approximants of PH 2 . 
Since the chain pattern is invariant with respect to space reflexions, it is sufficient to 
require that the solutions are bounded over a quarter of the chain. Rather than by 
diagonalizing products of exponentially large numbers of matrices, we have studied the 
spectral problem in the space of the traces, as this approach allows for a much higher 
precision. The boundedness condition quoted above translates into: 

i ■ \Trace{fy^ ■ f^y)\ < +1, ^ ■ \Trace{fy^)\ < +1. (12) 

The relevant property to be stressed here is that the space of the traces is hnite- 
dimensional. This is a direct consequence of the Cayley-Hamilton (CH) relation 

— Trace{M) ■ M + det{M)l = 0. (13) 

In fact it can be verihed that (CH) allows to reduce the order of positive powers of 
matrices M^, [k = 2,3,...) and to turn negative powers M~^ into positive ones. More 
importantly one can show that the traceless parts Ai of the matrices M 

M. = M — Trace{M) ■ 1 (14) 

satisfy Clifford anticommutation rules. This point can be easily ascertained by writing 
the M’s in the real representation 

M = mo • To + m • r, (to = 1, Ti = i ■ Us, Ts = as, = Ui). (15) 

This maps the matrices into a four-dimensional space, their coefficients being 
1 

rriQ =-Trace{M), m*, (i = 1, 2, 3). (16) 

The algebra of the spin operators r^, {i = 1,2,3) allows to compute, for any couple 
of matrices A, B, the anticommutator [H, R]_|_. By writing the resnlt in terms of the 
traceless parts A and B one hnally gets: 

A-B + B- A={A-B)-1, (17) 

where (H; B) = Trace{A ■ B) — \ - Trace{A) ■ Trace{B). 

These rnles allow to collect eqnal factors within prodncts, so that their powers can be 
reduced to order one. In this way we get monomials were the fonr types of vertices are 
raised at the most to power one. The occurrence of their T-inverted partners introduces 
in the algebra the matrix ai. We are thus led to operate in the larger space SL^{2, R). 
In the calculations we renamed the vertex matrices as follows: 

X(l) = Tyx, W(I) = Txy, -^(2) = Tyy, X (2) = Txx- (18) 

The elements W(l), X(I) correspond to the orthogonal vertices; prodncts where these 
factors occnr in alternating order are associated with ladders aligned with the diagonals 
of the sqnare lattice, as shown in Fig.4. In general the transfer matrix of a chain with 
n vertices has the form X(? 7 „, ..., 771 ) = X(? 7 „) •.. ■ X{rji) ■.. • X{rji) where the labels 
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rji run over the four values 1, 1, 2, 2. We wrote the trace map in terms of the following 
variables: 


x{ri,ri’,....) =Trace{X{ri,ri’,..)), (19) 

y{ri,r]’,...) = Trace(X(7],7 ]',■ Ui). 


3. Multifractal spectrum 

We will now illustrate the results obtained with the trace map for the scalar disordered 
laplacian A^;. We considered a very small anisotropy: ki = 1 — e, kt = (kl << !)• 
When e > 0 the x-variable at the horizontal links undergoes a restoring force weaker 
than at the vertical ones. In other words its longitudinal displacements are slightly 
softer than the transversal ones. 

One hnds that after one iteration the initial band is split into four bands. This branching 
process repeats itself at the successive steps so that the number of bands A4 scales as 
the size of the system A4 ~ 4^. 

We have found that upon increasing the anisotropy the bands start to overlap. The 
results reported here do not cover this case, as we made sure that the behavior A4 ~ 4^ 
was preserved up to /c = 10 iterations. 

The spectra of ’disordered laplacians’, when characterized by means of the density of 
states (DOS) p(A), can display a variety of scaling behaviors p(|A|) ~ (|A|)"“b Here 
instead of a we use the exponent /i that gives the scaling of the bandwidths AA with 
respect to Af: A A ~ Af~^. As illustrated in the caption of Fig.l p and a are simply 
related: fi = K In the literature on phononic spectra it is more usual to consider, 
rather than a or p, the spectral dimension dg [25] that gives the scaling of the density 
of states with respect to the frequencies p{u) ~ ( 0 ;)'^®“^. Since A = the two 
exponents are trivially related: a = y. Quite surprisingly we have found that close 
to the isotropic limit the spectrum has a whole set of power law singularities. In other 
words the spectrum is strongly sensitive to disorder even with rather small values of 
e. We did not sistematically explore the isotropic limit e —)■ 0 but the whole set of 
singularities [26] found here is qualitatively different from what would be expected with 
periodic chains. Here together with genuine normal modes, herefrom named Peano 
modes, corresponding to /i = 1, we hnd singularities with powers below and above 
/i = 1. Working on discretized, disordered extensions of laplacians such as the one 
examined here, even within a small spectral range AA a variety of singularities can be 
found, so that the numerical plots of the DOS versus A are generally hard to interpret. 
This difficulty can be surmounted when the singularities are of power-law type. In 
that case by collecting together the portions of the spectrum that share the scaling 
exponent one can obtain a smooth distribution. This procedure goes under the name of 
thermodynamic formalism [27] for multifractals. 

One dehnes a partition function Z{t), where r is the analog of an inverse temperature; 
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Figure 1. Isotropic limit of the disordered laplacian. We display here the distribution 
S{p) of the scaling exponents fj, obtained from the periodic approximant of order k=10. 
The three plots correspond to increasingly small values of the anisotropy parameter 
(e = 10“® red (circles), 10“® green (triangles), 10“^ blue (squares). It appears that 
the isotropic limit is characterized by a spectrum which turns out to be qualitatively 
different from what would be expected in the case of periodic chains (see text). We 
consider here the distribution for the exponent which defines the scaling of the 
bandwidths (AA « We recall that the more frequently used exponent a 

refers instead to the scaling of the density of states (DOS) p{X) « A““^. The two 
representations are equivalent; in fact from the boundedness of the integrated DOS 
H{X) : AH{X) = p{X) ■ AA = 0(1) one obtains p = 


the sum is taken over the J\f bands of a periodic approximant, and accounts for their 
individual widths AXi, (/ = 1, 

Z{t) = ^ (AA)'?W. (20) 

The exponent g(r) is the analog of a free energy: in fact it is the sum of an entropy 
and an energy contribution. To see this, let us consider the subset of bands E{fi) whose 
widths scale as AA ~ If we make the dependence on /r explicit on the r.h.s. of the 
equation, we get: 

g(r) = 5(/r) + /i ■ r, (/i = ^). (21) 

Fig.l is the plot of the ’entropy’ function S{fi) that weights the different /x-scalings 
within the spectrum, so that its value corresponds to the fractal dimension of the subset 
The extremal values ^iraax,min = hmT-^±oo ^ of the support of S{fi) characterize 
the ’vacuum’ (temperature T —0+) and the ’fully inverted state’ (temperature T —)■ 0_) 
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Figure 2. Spectrum corresponding to the vector-valued eigenstates: plot of the 
distribution S{^) obtained with inhomogeneity parameter e = 10“^ (green, triangles). 
For comparison we plot the distribution associated with the corresponding single 
variable problem (blue,circles). The interval {umim l^^max) appears to be contracted 
in the two-dimensional problem. Thus the requirement that the eigenvalues are at 
the intersection of the two single-variable spectra Xxikt, ki) and Xy(kt, ki) modifies the 
extremal regions ^ « fimin, M ~ fJ-max- 


regimes. Notice that from the formula relating q{T) with S{fi) one identihes the ’internal 
energy’ U of the system: 

dr 

So far we have considered the scalar operator A^., acting on the single variable x: this 
problem describes a ’polarized’ D=1 diffusion, as if a blind walker experienced a sequence 
of varying mobilities encoding the path in its meandering along the plane. The results 
obtained for the full vector-valued problem A(^ph) are shown in Fig.2, where they are 
also compared with the scalar case. The two distributions turn out to be different, in 
particular the range of exponents /i is larger in the scalar case. When the operator 
A(pjy) is diagonal in x and y, its eigenstates can be obtained by combining the solutions 
of the single variable problems in x and y: 

aA(n) = (a;A,(n),0)^ /3 aH = (0, |/A,(n))^, (23) 

where and Xy are the eigenvalues of A^, = Ax{kt, ki) and A^^ = Ay{kt, ki). In order 
to have an eigenstate of A(pp-) with eigenvalue A it must be: 



A Xyi^kfy ki^ Xxikt) ki'j. 


(24) 
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In other words since the spectrnm of ^[ph) is at the intersection of the two spectra it 
is not surprising to hnd a smaller range of exponents. 

4. Langevin equation 

We proceed now to examine the gaussian relaxation process associated with the 
disordered Laplacian A(^ph). This allows to understand the nature of the physical states. 
It turns out that the behavior of the fluctuations is determined by the spectral exponents 
/i. Our discussion relies on the literature dedicated to dispersion relations of fractals [28] 
as applied e.g. to the Generalized Gaussian Structures [29]. We assume that the vector 

held of deviations X(t) = .,a;(4^,f)) evolves according to 

the following Langevin equation: 

j^X{t)=V-A^PH)X{t)+m (25) 

< >=2-V ■ kpT ■ 6ij 6n,n', (hi = x,y). 

This process can be considered as a generalization of the Rouse model, describing the 
relaxation of an homogeneous polymer. The critical dynamics of self-avoiding Rouse 
chains and membranes has been studied by K.J.Wiese [30]. Here we are studying a 
purely linear system, but with a nontrivial dependence on the geometry of the underlying 
pattern. For the sake of simplicity let us omit from now on the vector indexes and put 
D = 1 and kpT = 1. It is well known that the process dehned by a Langevin equation is 
fully described by the response function (also called dynamic susceptibility) and by the 
correlation function; these two functions are related through the huctuation-dissipation 
theorem [31]. Before examining the case of the disordered laplacian it is useful to remind 
some basics on gaussian processes. We consider the relaxation of a held (j){x, t) in the 
presence of white noise ^{x,t): 

^0(a;, t) = -[r- V^]0(a;, t) + ^{x, t). (26) 

In the energy-momentum representation the dynamic susceptibility is x(A, Q) = {—i ■ 

A + r + Q2)-i. 

The correlation function is dehned as follows C{x,t) =< (f>{x,t) ■ 0(0,0) >= 

f ^ f j^g energy-momentum form is: 

G(A,Q) = 2-(A2 + (r + Q2)2)-i_ (27) 

We want now to relate the dihusive dispersion in the scale-free regime r = 0 with the 
band-edge singularity in the spectrum of the laplacian: more precisely we show that the 
dispersion is the counterpart of the multiplicity of states at the edges of the spectrum. 
In order to make this clear we consider the discretized version of the laplacian, where 
the counting of the states is elementary: 

^(0n+l + 0n-l - 2 • 0„) = Etpn, 


(0(iV) = 0(0)) 


(28) 
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Consider now the map E ^ n associated with the set of eigenvalnes E = E{n) = 
E{kn) = t{cos{kn) — 1), {kn = ^): this map is invertible with the exclusion of the band 
edges E = 0 and E = —2f, where the DOS displays the D=1 van Hove singularities 


dp{E’) 


1 dE' 

27r (E -E'^)^' 


(E' 


E + t)]. 


(29) 


The points n{E), when projected over the E-axis, produce in the interior of the band 
the scaling \dE\ oc \dn\ and \dE\ oc |(inp at the edges. Since one has \dn\ ~ Q, the 
latter scaling is the laplacian’s dispersion Q^. 

We generalize this connection to the case of the disordered laplacian: it is understood 
that our argument is merely heuristic because it refers only to the behavior in the scale- 
free regime, that is to the region A ~ 0 and Q oc N being the chain length. Recalling 
that within the spectral set E{p,) the singularity is |AA| ~ we associate to it a 

dispersion relation of the form 




(30) 


This amounts to saying that the contribution arising from E{fi) is described by an 
anomalous gaussian process undergone by a held 0^(A, Q) with dynamical susceptibility 

XM(A,Q) = H-A + |Qr)-^ (31) 


The correlation function 

2 

= (A2 + (|Q|/^)2) 

when integrated over A and Q, is represented as follows: 

C^{x,t) =< (pf,{x,t) ■ 0^(0, 0) >oc J (33) 

The infrared behavior associated with E{fi) is thus: 

C^(n, 0) ^ n » 1, (n(L) ^ L^), C(0,t) (t)^~^, t » 1. (34) 

Here we replaced the space variable x with the chain’s parameter n and L indicates the 
linear size of the planar region covered by n subsequent links. 

In conclusion, having shown that the exponent /i determines both space and time 
behaviors in the infrared regime, we can classify the different types of huctuations 
available to the system. 

The region p < 1 describes huctuations that decay in the length n as a power law. 
The corresponding eigenfunctions presumably have an oscillating behavior where the 
amplitudes are progressively reduced with distance, as in destructive interference. In 
addition, when p < 1 the huctuations are short-lived since they decay with time as a 
power law. 

In the region /i > 1 we hnd the complementary behavior, where the huctuations grow 
with n. They can be characterized in terms of the fractal dimension Df = or, in 
the polymeric interpretation, by means of the Flory exponent u = ^ = {p — l)/2. 

One can wonder how spatially growing functions could be acceptable eigenfunctions: the 
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point is that in the presence of oscillations the boundedness condition can be recurrently 
fulhlled. In this case there is presumably some sort of constructive interference. The 
two regions fi < 1 and fi > 1 are complementary also with respect to time dependence: 
in fact when fi > 1 the chain sites undergo subdiffusion with an exponent = 1 — K 
This behavior can as well be characterized by the dynamic exponent = 2/'Jd. The 
scenario found here leads to associate /i > 1 and p. < 1 with excitations respectively 
dominated by shear and by dilatation. 

It is worth mentioning that the fluctuation-dissipation theorem applies to the anomalous 
gaussian process described here. As a consequence the well-known deGennes relation 
= 2 + Df following from that theorem is fulhlled for all exponents /i. The case 
fi = 1, which marks the separation of the two regions described above, corresponds to 
oscillations with stable amplitudes. The corresponding Peano modes will be discussed 
in connection with the symmetries of the system. 

We conclude this Section with a comment on the general case of chains PH a. Being 
always generated by deterministic inhation, it is reasonable to expect that their self¬ 
similarity should be mirrored into nontrivial spectral properties. The eigenvalue problem 
in D dimensions, if approached along the lines presented in this paper, would involve 
transfer matrices of dimension (2D x 2D). 

This can be difficult to deal with even in the case D = 2: we remind that it is possible 
to improve the numerical resolution, by using traces instead of matrices, only when the 
components x and y are decoupled. One can try to overcome this complexity by relying 
on a mean-held approach, which in the present context amounts to disregarding the 
huctuations of fi. Taking into account that the chain is embedded in dimension D, self- 
consistency suggests that the mean-held exponent fi(D) should be chosen in such a way 
as to make the fractal dimension Df to coincide with the bulk’s dimension; Df = D. 
In order to determine fi(D) let us consider a couple of points (1,2) placed at a distance 
r in the bulk, and use the formula Df = In this way we hnd: 

< (T(l) — T(2))^ >cs = no = l^{D) ^ ^ 

Notice that in the language of membranes the huctuating held x(n, t) undergoes the 
roughening transition at /i = 1 , hence the mean held sets the membrane in the rough 
phase [32]. 

In the planar case we obtain fi(2) = 2 = Df. It turns out that huctuations of exponent 
fi = 2 would correspond to structures softer than the elastic network of PH 2 itself. 
In fact the maximal growth of the latter is determined by fimax and we have found: 
^^max < h( 2 ) = 2 , (7 = 0 , e —)■ 0 ). 

In D=3 the mean held exponent is /i(3) = |, which implies the anomalous dihusion 
C(0,t) ~ P. This behavior was actually found to occur with high probability in 
experiments on bacterial chromosomes [7, 33] and is compatible with conhgurations 
close to the optimal packing Df = 3[9]. Let us hnally notice that huctuations having a 
fractal dimension higher than D cannot be excluded. In fact in some of the measurements 
quoted above [7] also rather small anomalous dihusion exponents 7 ^ ('yo < f) were 
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dq/dx 



Figure 3. Plot of the function ^(r) obtained from the bands of the approximant 
of order k = 10. The inhomogeneity parameter is e = 10“'^. The ’internal energy’ 
U of the system is given by U = — The asymptotic values of the function are 
respectively Umim at large negative values of r, (r —>■ —oo), and fimax, (r —>• +oo) at 
large positive values of t. At the critical (second order) point e —>■ 0 the function ^(t) 
has a single inflexion at r = 0. The plot displayed here clearly shows three inflection 
points, indicating the onset of the phase separation: the ’internal energy’ U is going 
to develop two minima. 


obtained: they would imply 1 < p < /i(3) and Df > 3. 


5. Inverse energy cascade and phase separation 

In this Section we describe the global properties of the system’s relaxation by relying 
on the results discussed above. We have found that the thermodynamic formalism of 
multifractals, which up to now has been used primarily for technical reasons, involves 
a physical interpretation of the ’internal energy’ U = — ^(r) that is consistent with 
the scenario of relaxation. The region of negative r describes the regime of negative 
temperatures (i.e. the inverted system) and in fact, as one can see from Fig.3, there one 
finds the higher values of U. The limit r —)■ +cxo corresponds to the ground state of the 
thermodynamics and there the ’internal energy’ U has its minimum. 

The interesting point, revealing some connection of U with the energies A, is that the 
infrared region A ~ 0 appears to be dominated by the largest scaling exponent fimax > 1: 
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at least this is what emerged from our numerical checks. As discussed in the previous 
Section the values /x > 1 are associated with excitations of a large size, presumably 
dominated by shear. In the opposite case, i.e. whenever the dilatation dominates, one 
expects localized excitations and higher energies, associated with values /i < 1. 

If the ’internal energy’ U is monotonic with respect to the energy content of the 
fluctuations, then from the plot reported in Fig.3 one gets a clear picture of the process 
of relaxation: when prepared on localized states, the system is led to evolve towards 
the shear-dominated regime /i > 1. This behavior has been in fact observed in two- 
dimensional colloids [17, 18] and is intriguingly similar to the inverse energy cascade of 
two-dimensional turbulence [ 22 ]. 

It turns out that the analogy between this model and two-dimensional colloids is not 
limited to the space-time behavior of the excitations, but involves also the spectral 
properties. In fact if one were to add to the model the contact couplings among adjacent 
but nonsubsequent sites, the new network would describe a sort of textured membrane. 
The resulting phononic spectrum, in addition to the contributions arising from the Peano 
pattern, would display in the density of states p(A) a contribution p(A) d =2 ~ A, arising 
from planar phonons, as in the standard Debye spectra of two-dimensional systems. The 
quoted experiments on colloids did in fact display at low frequencies an anomalously 
high density of states, superimposed to the D=2 Debye behavior. 

So far we have discussed the spectral properties in the closeness of the isotropic limit, 
where the spectral bands are predominantly separated. The question then arises about 
how the scenario changes when the anisotropy is large enough to make the bands overlap. 
We have found that the anisotropy induces a phase transition [34]: the spectral measure 
S'(/i) progressively splits into two separate distributions. Possibly the simplest example 
of a complete phase separation in a spectral measure is provided by the D=1 periodic 
chain. As we earlier recalled, in that case there are only two exponents: n = 1 and 
/i = 2 . Since their weights are S'(/i = 1) = 1 and S'(/i = 2) = 0, the periodic chain 
is located at the extremal point of a line of phase coexistence, at the extinction of the 
phase = 2 , where the system is fully converted into the jj, = 1 phase. 

A phase coexistence was fonnd [35] in the spectrum of the Quantum Ising Quasicrystal 
(QIQX). The QIQX model describes a chain of spins where two Ising couplings are 
aligned in quasiperiodic order [36]. When the two conplings are nearly equal (Jo ~ Ji) 
the function S{fi) has snpport /imin < < Umax with ^min = 1 and ^max = 2. In other 

words the spectrum is at the coalescence of the two mentioned phases of the periodic 
chain. As the two conplings are made progressively different (Jq ^ Ji) the two phases 
separate. 

In general a phase separation can be hardly detected by inspecting S'(/i), because the 
numerics naturally converge to the convex envelope of the two distributions. In fact 
S{fi) is the Legendre transform of the object of relevance in this context, which is the 
’internal energy’ \] {U = By closely inspecting Fig.3, where we display ^ at 

the value e = 10 “^ of the anisotropy parameter, one can notice three inflection points: 
A ^ (—4,.!.), B ^ (0.,1.1), C ^ (2., 1.3). At the (second order) critical point one 
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Figure 4. Planar strip associated with products where the two types of orthogonal 
vertices -^(1), -^(1) are ordered in alternating fashion. The horizontal (vertical) steps 
of the ladder correspond to discrete displacements along the space coordinates qx, (qy)- 
The ladder can be associated with a trajectory undergoing a sequence of scatterings at 
the strip’s boundaries: the corresponding motion, aligned with the direction depicted 
in black, would carry a linear momentum p in that direction. This suggests that the 
orthogonal vertices ^(1), -^(1), which represent the turning points of the chain, should 
be associated with displacements in the plane of momenta {px,Py) where the axes form 
an angle 7r/4 with respect to the orientation of {qx,qy). 


expects a single inflexion point, here we are very close to that point. Upon increasing 
the parameter e a local maximum is expected in the region delimited by {A,B). This 
maximum identihes the emergence on the left side of the plot, i.e. in the inverted regime, 
of a phase whose dominant scaling evolves towards pmm- On the right side the scaling 
of the other phase evolves towards ^max- It is not difficult to associate the former with 
the solid, where the high energy dilatation modes are dominant, and the latter, where 
the fluctuations grow both in n and in t, with the shear-dominated liquid. 


6. Critical crystal with eightfold symmetry 

In the previous Section we have shown that the model undergoes a phase transition, and 
that its isotropic limit is found at the saddle separating two phases, respectively of solid 
and liquid type. The Peano network can thus be considered a critical crystal, marking 
the point where phase coexistence turns into coalescence. This situation appears to be 
the D=2 analog of the well studied critical regime at the D=1 delocalization transition 
[23, 37]. It is thus appropriate to remind some results regarding that case. 

The D=1 transition has been described by Wilkinson within a semiclassical (WKB) 
approach based on the universal hamiltonian [38]: 

-■ H = cos{p) + a ■ cos{q) (36) 

Exactly at the critical point a = 1 the hamiltonian has a twofold symmetry, since it 
is self-dual in the exchange of the (adimensional) space and momentum variables. It 
identihes a square lattice in the (g,p) plane, the lattice points being at the saddles where 
the separatrices meet. These separatrices are soliton trajectories of energy 77 = 0; apart 
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from them, every other orbit is closed, so that when a = 1 the particles are conhned 
both in the q and p directions. 

When the self-duality is broken, e.g. when a > 1, aligned solitons merge into trajectories 
running parallel to the p-axis: in this region the particles are localized with respect to 
the g-direction. In correspondence with the complementary coupling (^) one hnds 
the dual trajectories, where the particles are instead extended in the g-direction. 
It resulted from Wilkinson’s analysis that the critical properties obtained with the 
universal hamiltonian, which has a twofold symmetry, actually hold for an entire class 
of hamiltonians, characterized by fourfold symmetry. The following is an example of 
such hamiltonians: 

1 

-H = cos{q) + cos{p) + cos{2 ■ q — p) + cos{q + 2 ■ p) (37) 

We suggest that this fourfold symmetry could originate from the operators representing 
the discrete translations in the g and p directions, together with their time-reversed 
counterparts. Products of these four elements are identihed by the oriented paths on the 
lattice phase space. In order to see how the D=1 case is related with our construction it is 
useful to examine the quantum hamiltonians based on chains organized according to the 
two-letter (A, B) Fibonacci substitution rule. The Fibonacci chain is often depicted as a 
ladder over a two-dimensional lattice, where the horizontal and vertical steps correspond 
to the letters A and B. The steps are combined in such a way as to £t the ladder within 
a strip having slope m = ■ The ladder shown in Fig.4 is represented by products 

where the two operators W(l) and W(l), associated with the orthogonal vertices of our 
construction, are ordered in alternating fashion. In this elementary configuration the 
ladder is contained within a strip of slope m = 1. One can notice that since the operators 
X(l) and W(I) involve couples of steps {AB) or {BA), they are more properly associated 
with translations in the momenta along the diagonal. The quantum mechanics on the 
ladder has the form of a constrained D=2 hamiltonian system where translations in the 
momenta are reduced to a single dimension. 

The fact that at the critical point the D=1 quantum system is correctly described by 
the WKB approximation translates, in the language of trace maps, into the condition 
Trace{A ■ B ■ A~^ ■ B~^) = 2, meaning that the group commutator is the identity. This 
condition is connected to the localization properties of the wave functions through the 
following trace relation: 

Trace{A ■ B ■ A~^ ■ B~^) = — xyz — 2 (38) 

X = Trace{A), y = Trace{B), z = Trace{A ■ B). 

On the r.h.s of the above one recognizes the well-known constant of motion of the 
trace map [19] representing the surface I{x,y,z) = x"^ + y"^ + z"^ — xyz. The topology 
of this surface changes when J = 4, i.e. when the group commutator is trivial, and 
correspondingly two separate types of map trajectories (localized-delocalized) meet. 

We now return to the trace map of our case, leaving to future work [39] a 
detailed analysis. The map involves the monomials made out of the four vertices 
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Figure 5. Graphical representation of the elementary operations acting on the four¬ 
dimensional lattice. Upper line: the blue (red) square path represents the loop of 
discrete translations in the p-plane (q-plane). Lower line, right side: the eightfold 
motif obtained by projecting the p-loop on the plane of the q-loop. 

Lower line, left side: the loop formed by combining the eight operators that generate 
the Peano chain. The colour convention is consistent with the rest of the Figure. The 
red lines represent paths in the q-plane, hence the red perimeter describes a loop in 
the q-plane. The black dots represent the vertices: the flat vertices ^(2), X{2) carry 
only red lines, the orthogonal vertices ^(1), -^(1) are decorated with blue lines. The 
blue segments are aligned with the translations on the p-lattice, as illustrated in Fig.4. 
If the blue segments are continued up to their intersections, the eightfold motif shown 
on the right side is reproduced. 


(X(l), ^(1), -^(2), X{2)). One has six order-two monomials, four order-three 
monomials and one order-four monomial. The degrees of freedom are doubled since 
the four basic vertices occur with their T-reversed copies. The map is thus dehned 
over the traces x = Trace{X) and y = Trace{X ■ ai) illustrated in Section 2. The 
overall space has dimension 15X2 but this large number can be reduced by exploiting 
some symmetries. For instance the invariance under cyclic permutations generates a 
triangular symmetry for cubic monomials. Cubic monomials are further related through 
some trace identities[39]. 

On top of that the occurrence of 4 x 2 basic vertices enforces an eightfold symmetry. 
This, we believe, appears as the extension of the D=1 case that is characterized by a 
fourfold symmetry. In fact in quantizing an hamiltonian system over a (2 x 2) phase 
space one needs (4 x 2) operators representing the discrete {qx,(ly), (Px-iPy) translations 
and their reversals. These operators act in this case on the vector of elastic displacements 
having values in a real two-dimensional vector space. 

In order to understand how this geometry effects the conhgurations of the system, we 
project the lattice of momenta on the plane of conhgurations. In Fig.5 we enumerate 
the discrete translations of the phase space and compare them with the algebra of 
the vertices. The motif displayed on the lower left is obtained by combining the 
vertices X(l), X(i),X(2), X(2) with their reversed counterparts. On the lower right 
we show the effect of projecting the {px,Py) lattice over (g^,, qy). Notice that in the usual 
representation of the union of a q-lattice with its dual one assigns the nodes of the dual 
to the centers of the q-plaquettes. Under that convention the tiling of the plane would 
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Figure 6. The octagon and its descendants: this figure illustrates how the octagonal 
symmetry is compatible with dilatations. Discrete rotations and discrete dilatations 
match, as in regular crystals rotations match with lattice translations. This ’polar 
crystal’ property is shared by the Peano chain with the octagonal quasicrystal. 
Periodic approximants of the octagonal quasicrystal are obtained by projecting the 4D 
hypercube on the plane; the resulting tiling turns out to have an average connectivity 
substantially higher than in the Peano network. 


appear in the form of two equally oriented square lattices. The eightfold motif obtained 
here results instead from the convention of making the centers of the q and p plaquettes 
to coincide and to having them rotated of an angle 7r/4. 

We hnd that this convention reproduces the geometry associated with the operators 
X(l), X(1),X(2), X{2). In fact the two motifs displayed on the second line give rise to 
identical periodic tilings of the plane. 

It turns out that in recent experiments on bilayer water, precursors of crystallization 
were observed to generate quadratic, hexagonal, pentagonal patterns, and more 
importantly, also the dodecagonal quasicrystal [24], On very general grounds, since 
quasicrystals are much softer than solids, it is somewhat to be expected that in the 
liquid to solid transition various structures of intermediate stability should be explored 
by the system. 

Given this eightfold symmetry the question is then: to what extent is the Peano network 
related with the octagonal quasicrystal? The effective coordination number of the 
network is certainly smaller than the average coordination number of the octagonal 
quasicrystal. Periodic approximants of the latter [40], obtained by projecting the 4D 
hypercube on the plane, have square-shaped elementary cells. Each cell has vertices 
with coordination numbers z varying from three to eight {z = 3,4, 5, 6, 7, 8). The 
free particle (laplacian) quantum motion dehned over such tilings is characterized by 
anomalous diffusion, but with exponents much larger than those obtained here [41]. The 
reason for that can be easily understood, since in the octagonal quasicrystal there is a 
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Figure 7. Peano modes at the matching conditions of the ’polar crystal’. Plot of the 
number of solutions M versus the energy A in the closeness of A = 0; the distribution 
is defined over intervals having width 6X = 10“®. Apart from the normalization, the 
plot is a coarse-graining of the density of states (DOS) p(A). The data refer to the 
approximant of order k=6 and to the anisotropy parameter e = 10“®: we are thus quite 
close to the critical point. The points fall along five curves: Af{l, A), (Z = 1, 2, 3,4, 5). 
The curves a = cr(A, s) orthogonal to the AA’s are coloured in red. The a curve 
intercepts approximants of increasingly large orders 1. 


larger number of channels open to hopping. 

We would like to stress that the octagonal symmetry arising in the Peano network 
does not originate in the euclidean metrics, but in the hamiltonian structure. If the 
system is integrable the hamiltonian can be made quadratic and in the isotropic limit 
the rotational invariance implies the following form: 

H = a-{pl+pl)+(3-{ql + ql). (39) 

At the delocalization transition the system is self-dual, thus at the critical point it must 
he a = (3\ 


ipl + pD = (9^ + i) 


(40) 
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A 




Figure 8. Lower order symmetries compatible with the octagon. One can identify 
the fourfold loop in the p-plane (in blue) A'B'C'D', and, in red, its analog in the 
q-plane ABCD. Triangular loops mix the q and p coordinates. We show here two such 
loops: AB'C and CD'A'. On the right side we display two fivefold loops obtained by 
connecting the vertices left out of the triangular ones. The loop depicted in green is the 
complement of AB'C. Within a semiclassical description the loops can be associated 
with scattering processes adding up to zero momentum exchange, as required by scale 
invariance at the critical point. 


The above means that the allowed states are points of the (2 x 2)-dimensional lattice 
touched by the sphere of the energies and satisfying the self-duality condition. 

This allows to understand why the Peano chain can sustain normal modes, in spite of 
its intricate structure. As shown in Fig.6 the octagon reproduces itself at discrete scales 
where the dilatations match with discrete rotations. The chain shares with the octagon 
this ’polar crystal’ property and can thus support normal (Peano) modes, organized in a 
hyerarchical fashion. The numerical results conhrm this argument. In fact by inspecting 
the multiplicity of solutions as a function of the energies A displayed in Fig.7 one can 
hnd some relevant features enumerated below. 

A)The distribution appears as the union of several curves 7\/'(/, A), (/ = 1,2, 3,4, 5) 
having shapes similar to what one would obtain with periodic chains. The curves can 
be associated with approximants of increasing order 1. In fact they represent solutions 
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of growing multiplicity, as one can see from their intercepts with the vertical axis. 

B)A second family of curves a (in red) is apparently orthogonal to the J\f{l,X). 
The (t’s can be identihed by means of their intercept A with the horizontal axis: 
a = a{\,s), \\ —)■ (t(A,s = 0 ). As the arc parameter s grows the curve o-(A,s) crosses 
approximants of increasing order /: the multiplicity (degeneracy) of the eigenvalue A 
increases accordingly. 

The Peano waves can be depicted as Rayleigh waves whose packet is split into non- 
dispersive contributions corresponding to the different scales of approximants. This 
scale separation actually involves not only the normal modes (/i = 1 ), but also the 
anomalous modes 7 ^ 1: one can identify the /—curves down to the edge A ~ 0. 

In addition to the octagonal symmetry, lower order ones, enumerated in Fig. 8 , can as 
well be of importance. The case of triangular symmetry was shortly illustrated above, 
in association with cubic monomials. 

7. Conclusions 

We have studied the elastic deformations of the plane-hlling closed Peano chain, taken 
as an elastic network. The shear and dilatation modes are intrinsically coupled, because 
their waves are constrained to propagate along a thin planar strip undergoing the se¬ 
quence of turns dictated by the deterministic pattern of the chain. 

The situation is thus similar to the interferences at the surface of an elastic material, 
giving rise to the dispersive Rayleigh waves. The Peano chain has the property of being 
self-similar, with periodic approximants. In a situation of this type the Rayleigh waves 
can hnd resonance conditions necessary for coherent propagation. The corresponding 
normal modes, named Peano modes, occur when the rescaling associated with the in¬ 
flation procedure matches with the discrete rotational symmetry of the pattern. 

We have studied the properties of the spectrum by means of the so-called thermo¬ 
dynamic formalism for multifractals. In addition to the Peano modes we have found 
solutions that can be qualitatively described as oscillating functions modulated by am¬ 
plitudes which grow or decay: they are associated with gaussian relaxation processes 
characterized by a scaling exponent /x 7 ^ 1 . 

These anomalous modes result from the frustrated coupling of shear and dilatation: 
when the resonance conditions are not exactly fulhlled, the waves can be amplihed or 
damped, depending on how interference operates. When /i < 1 the equal time correla¬ 
tion function decays along the chain as n^~^, hence these states are of localized type; 
correspondingly the chain sites undergo a damping process, with the power law t 
The fluctuations characterized by /i > 1 have a complementary behavior: rather than 
damping they undergo anomalous diffusion and their equal time correlation functions 
grow in amplitude along the chain. They are predominantly found at low energies 
(A ~ 0). Taking the above into account, the relaxation of the system can be described 
as a flow from localized high energy structures towards extended low energy ones, where 
the latter survive in the long time regime of the process. This picture turns out to repro- 
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duce, at least qualitatively, the behavior observed in two-dimensional colloids [17, 18]. 
It may then be not coincidental that also the presence, observed in the same exper¬ 
iments, of a low-energy bosonic peak can be readily explained within our model. In 
our opinion the peak reveals a contribution to elastic energy arising from closed, large¬ 
sized one-dimensional patterns. Structures of this sort are not easily described within 
a membrane-like network model. The density of states should then account for these 
excitations, which add up to the standard D=2 phononic part. 

Possibly the main result obtained in this work is the proof of the existence of a critical 
point in the isotropic limit of the model, where two phases coalesce. These phases are 
respectively dominated by the /i < 1 and /x > 1 regimes: localized versus delocalized, or, 
in other words, solid versus liquid. It is thus appropriate to consider the Peano network 
as a critical crystal. 

Since, as previously mentioned, the dilatation modes have higher energies and are more 
localized, the energy flows towards extended configurations and lower energies, as in 
the inverse energy cascade of two-dimensional turbulence. We have argued in addition 
that this critical point should be characterized by an eightfold symmetry: from this 
perspective the Peano chain is closely related with the octagonal quasicrystal, but has 
a much lower connectivity. Quite intriguingly, experiments on bilayer water revealed a 
transition from the liquid to the dodecagonal quasicrystal[24]. The eightfold symmme- 
try that emerges here is strictly related with us limiting the description to contributions 
from large-sized D=1 structures; enlarging this picture could lead to results of a greater 
relevance for the phenomenology, but at the expense of losing the solvability. 

From the perspective of the possible connections with the spatial ordering of chromo¬ 
somes, the mean-held approach presented in Section 4 gives an estimate of the anomalous 
diffusion exponent that hts quite well with the numbers observed in bacteria. 

Let us hnally comment on low energy excitations of a different kind that should as well be 
expected in Peano-shaped structures. The chain phonons treated in this work describe 
the small deviations with respect to a fully ordered pattern, but in addition one could 
consider isolated topological defects, giving rise to deviations from exact space-hlling. 
The motion of such defects necessarily involves collective reorderings along the chain, 
i.e. phasons [42]. We are led to associate with excitations of this kind the correlated 
motions observed in chromatin [43]. At even higher energies, multiple defects should be 
accounted for, but this goes far beyond the perspective of our approach; nonetheless, 
we point out that the chain phonons should give a relevant contribution in the infrared 
regime, since they carry the smaller momenta available to the system, scaling as L~^. 
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Figure Al. The convention used for labeling the different types of vertices is here 
specified. As the iteration is based on a quadripartite structure, in order to assign 
labels to vertices one can start with the elementary square-shaped closed path. Since 
paths can as well go through lattice sites without bending, in addition to orthogonal 
vertices also the flat ones (horizontal and vertical) must be accounted for. With the 
numerals (1,2,3,4) and (5,6) we respectively classify the patterns of orthogonal and 
flat vertices in the clockwise path; the vertices corresponding to the anticlockwise path 
are labeled with (1,2,3,4) and (5,6). 


Appendix 

In this Appendix we determine the symbolic substitution rule allowing to generate 
approximants of the closed Peano-Hilbert chain PH 2 . To the best of our knowledge 
the procedure illustrated here is original. It is organized in such a way as to preserve 
the correct matching of the four quarters making up a closed chain, as required in the 
spectral problem examined in the paper. The approximant of order k of PH 2 is a 
path of length Nk = 4^ that connects with no intersections all points of the (2^ x 2^) 
square lattice. We display in Fig.A2 one of the mentioned quarters, corresponding to 
the approximant of order k=4. Since the plane has the topology of a closed path 
cuts the plane in two polar regions. At first sight it is equivalent to assign a model over 
the square lattice or over a Peano approximant, but the two choices are topologically 
different. In fact the lattice is a torus while the Peano chain has the topology of a circle 
(see also the caption of Fig.A3). 

Let us now enumerate all possible vertices, in order to identify chain patterns with 
symbolic words. We label the available vertices with the numerals n and h, (n = 
1,2, ..,6), where n and h are oppositely oriented as illustrated in Fig.Al. The 
approximant of order k is thus written as a symbolic word made out of the 6X2 
letters I of the alphabet dehned above: 

= l{Nk)l{Nk - 1)./(3)/(2)/(l). (A.l) 

The alphabet used for the vertices is exceedingly extended, as every vertex is the image 
of a flat or orthogonal one under suitable transformations. These transformations are 
the discrete rotation R, the inversions P^, Py with respect to the x and y axes and the 

path inversion T. They act as listed here: T(m) = m, T(m(n)....m(l)) = m(l). m{n)] 

-Pa;(4) = 3, ^^(1) = 4, Py{Px{2)) = 4; 4 = R{1), 1 = R{2). We thus single out as basis 
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Figure A2. Closed PH 2 chain, also known as Moore curve. The lower left quarter 
{L C) of the approximant of order k=4 is displayed here. The arrows depict the 
ingoing and outgoing links connecting this quarter to its neighbours. The connection 
to the neighbour lying on the right side is provided by an horizontal link placed at 
the uppermost available location. A complementary situation regards the connection 
with the upper neighbour: in this case the link is vertical and placed at the leftmost 
available location. The inflation procedure is defined in such a way as to guarantee 
that this structure is preserved at all orders. 




Figure A3. In the lower line of this figure are displayed the two choices available 
when inflating the vertex to = 4; the parent (inflated) patterns are depicted in blue 
(red). With the left choice, the new vertices formed at the entrance and at the exit of 
the TO = 4 motif are respectively flat horizontal (5) and orthogonal (3). The convention 
chosed in the paper is displayed on the right: an horizontal link enters the orthogonal 
vertex (1), and a vertical link exits the flat vertical vertex 6. In the upper part of 
the figure we display the two resulting closed patterns after one step of the inflation 
procedure; the initial configuration is the squared perimeter depicted in blue. The two 
complementary H-shaped portions of the plane enclosed by the red patterns, when 
glued together as the two patches of a tennis ball, form a compact surface 8“^. Hence 
the approximants separate the compactified plane in two polar regions. The Peano 
chain is the maximally extended ’equator’ of this partition. 


elements the vertices m = 4 and m = 6. Our aim is to write the substitution rule under 
the constraint of closed path. Quite remarkably the constraint is sufficient to identify 
the rule. Let us consider e.g. the choices available for inflating the vertex m = 4: as 
shown in Fig.A3 two choices are available and they would imply twofold branchings at 
each inflation step. This ambiguity is eliminated if one requires that the inflated path 
must be closed. Every approximant can thus occur in only two forms; they are displayed 
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in the upper line of Fig.A3 for the first approximant. In our calculations we have chosen 
the option shown on the right column of Fig.A3. With this choice the four quarters 
making up the chain are connected through the two external (left and right) vertical 
links and the two central horizontal links: the H-shaped pattern displayed on the right 
side of Fig.A3 reproduces itself at all orders. Going back to the inflation rule, we write 
the convention as follows: 4 4' = 4 (6431) and 6^6' = (4231). Considering the 

symmetries of the pattern, the following identities must hold at every order: 

1 = TPy(4), 3 = TR(4), 2 = R-^Py(4), 1 = rR-\4), 3 = TP.(4).(A.2) 

We can thus write the inflation in terms of the two words Q 4 and Q^: 

Qe'’ = Qf~^\R-^PyQf~^'>){TRQf~^'>){TPyQ^~^^) (A.3) 

= Qt^Qf\TP,Qt^){TR-^Qt^). 

The word determines the lower left quarter of the closed chain The 

latter is obtained by operating with the symmetries and has the form: = 

{TP^^) {P.PyQf ) {TPyQf.) 
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